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Abstract 

We have performed a dynamical coupled-channels analysis of available p(e, e'Tr)N data in the 
region of W < 1.6 GeV and Q 2 < 1.45 (GeV/c) 2 . The channels included are 7*iV, ttN, rjN, 
and ttttN which has ttA, pN, and <jN components. With the hadronic parameters of the model 
determined in our previous investigations of ttN — > ttN, ttttN reactions, we have found that the 
available data in the considered W < 1.6 GeV region can be fitted well by only adjusting the bare 
7*iV — > N* helicity amplitudes for the lowest N* states in P33, Pn, Su and D13 partial waves. The 
sensitivity of the resulting parameters to the amount of data included in the analysis is investigated. 
The importance of coupled-channels effect on the p(e, e'Tr)N cross sections is demonstrated. The 
meson cloud effect, as required by the unitarity conditions, on the 7*iV — > N* form factors are also 
examined. Necessary future developments, both experimentally and theoretically, are discussed. 

PACS numbers: 13.75. Gx, 13.60.Le, 14.20.Gk 
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I. INTRODUCTION 



The electromagnetic parameters characterizing the excited nucleons (N*), in particular 
the 7*iV — > N* form factors, are important information for understanding the hadron struc- 
ture within Quantum Chromodynamics (QCD). With the efforts in recent years, as reviewed 
in Ref. the world data of 7*iV — > A(1232) form factors are now considered along with the 
electromagnetic nucleon form factors as the benchmark data for developing hadron structure 
models and testing predictions from Lattice QCD calculations (LQCD). The main objective 
of this work is to explore the extent to which the available p(e, eV)iV data in W < 1.6 GeV 
can be used to extract the 7*iV — > iV* form factors for the iV* states up to the so-called 
"second" resonance region. 

We employed a dynamical coupled-channels model developed in Refs. 0, 0, 0, 0, 0] . This 
work is an extension of our analysis [4] of pion photoproduction reactions. We therefore will 
only recall equations which are relevant to the coupled-channels calculations of p(e, e'n)N 
cross sections. In the helicity-LS J mixed-representation where the initial 7iV state is specified 
by its helicities A 7 and A at and the final MB states by the (LS)J angular momentum 
variables, the reaction amplitude of 7*(<f, Q 2 ) + N(—q) — > ir(k) + N(—k) at invariant mass 
W and momentum transfer Q 2 = —q^q^ = q 2 — uj 2 can be written within a Hamiltonian 
formulation [2| as (suppress the isospin quantum numbers) 

T J LSM , XN {k, q, W, Q 2 ) = t J LSMjXN (k, q, W, Q 2 ) + t^ N7TNMXN (k, q, W, Q 2 ) , (1) 

where Sn = 1/2 is the nucleon spin, W = u + E^{q) is the invariant mass of the 7*iV 
system, and the non-resonant amplitude is 

(k,q,W,Q 2 ) 

= v ls N 77N,x^\ n (^^ Q 2 ) + ^2 ^2 k' 2 dk't'l SNnNL , SIMIBI (k,k',W) 

M'B' US' 

xG M/B ,(k',W)vl, s (k>,q,Q 2 ). (2) 

In the above equation, Gm'b'^'i W) are the meson-baryon propagators for the channels 
M'B' = tcN, rjN, it A, pN, aN. The matrix elements "/sMBA^iMiQ 2 )! which describe 
the 'jN — > MB transitions, are calculated from tree-diagrams of a set of phenomenological 
Lagrangians describing the interactions between 7, it, r], p, lu, a, N, and A(1232) fields. The 
details are given explicitly in Appendix F of Ref. [2| . The hadronic non-resonant amplitudes 
IlSnkN L'S'M'B'(k-> k' 1 W) are generated from the model constructed from analyzing the data 
of ttN — > tcN, tctcN reactions [1, 0] . 
The resonant amplitude in Eq. ([T|) is 

t^N^xA^ ?. w, Q 2 ) = E [r J N*,L SN « N (k, w)]*D id (w)r J N * XjXN ( q , w, q 2 ) , (3) 

N*,N* 

where the dressed iV* — > nN vertex Fff*^ LSNnN (k, W) and A^* propagator D it j(W) have been 
determined and given explicitly in Ref. j5J. The quantity relevant to our later discussions is 
the dressed 7*iV — > N* vertex function defined by 

r^A^^Q 2 ) = r^, A7AiV (g,Q 2 ) 

N*, US' M'B' 

(k',W)G M 'B'(k',W)vl, s (k',q, Q 2 ) . 

M'B' US' 
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FIG. 1: Graphical illustration of the contribution to the ttN intermediate state to the dressed 
7*iV —> N* vertex defined by Eq. (gj). 



The second term of Eq. (j4]) is due to the mechanism where the non-resonant electro- 
magnetic meson production takes place before the dressed N* states are formed. This is 
illustrated in Fig. [1] for the contribution due to the M'B' = ttN intermediate state. Similar 
to what was defined in Ref . 0, S] , we call this contribution the meson cloud effect to define 
precisely what will be presented in this paper. We emphasize here that the meson cloud 
term in Eq. is the necessary consequence of the unitarity conditions. How this term and 
the assumed bare N* states are interpreted is obviously model dependent. This issue as well 
as the questions concerning the extractions of form factors at resonance pole positions will 
be discussed elsewhere, and will not be addressed here. 

Within the one-photon exchange approximation, the differential cross sections of pion 
electroproduction can be written as 



da 5 



dE e >dtt e /d{l* 



a T + ea L + J2e(l + e)cr Lr cos0* 



+ea TT cos 20* + h e J 2e(l — e)a LT ' sin</>* 



(5) 



Here T 7 = [a/ (2n 2 Q 2 )}(E e ' / E e ) [|?l|/(1 — e)]; e is defined by the electron scattering angle 9 e 
and the photon 3-momentum qi in the laboratory frame as e = [l+2(|g^| 2 / Q 2 ) tan 2 (0 e /2)] _1 ; 
h e is the helicity of the incoming electron; 0* is the angle between the 7r-iV plane and the 
plane of the incoming and outgoing electrons. The quantities associated with the electrons 
are defined in the laboratory frame. On the other hand, structure functions of 7*iV — > ttN 
process, a a = a a (W,Q 2 , cos 9*) (a = T, L, LT,TT, LT'), are defined in the final ttN center 
of mass system. The formula for calculating a a from the amplitudes defined by Eqs. CQ)-© 
are given in Ref. [^]. 

In this first-stage investigation, we only consider the data of structure functions a a of 
p(e,e'ir°)p 0, EH and p(e, e'ir+)n Ej up to W = 1.6 GeV and Q 2 = 1.45 (GeV/c) 2 . 
The availability of the data in the corresponding (W, Q 2 ) region are found in Table [H The 
resulting parameters are then confirmed against the original five-fold differential cross section 
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TABLE I: Available structure function data at Q 2 < 1.45 (GeV/c) 2 . 
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a The data are available up to W = 1.51 GeV 
b The data are available up to W — 1.41 GeV 



data [l4J. This procedure could overestimate/underestimate the errors of our analysis, but 
is sufficient for the present exploratory investigation. 

In section II, we present the results from our analysis. Discussions on future developments 
are given in section III. 



II. ANALYSIS AND RESULTS 

To proceed, we need to define the bare 7*iV — > iV* vertex functions T^* A A 

MQ 2 ) of 

Eq. (Tj0). We parameterize these functions as 



where qn and go are defined by M^* = q^ + E^{qR) with N* mass and W = go + -EW(go)) 
respectively, and 

G X (N*,Q 2 ) = A X (N*,Q 2 ), for transverse photon, (7) 
= S\(N*,Q 2 ), for longitudinal photon. (8) 

For later discussions, we also cast the helicity amplitudes of the dressed vertex Eq. (j3j) into 
the form of Eq. (jSj) with dressed helicity amplitudes 

A X (N\Q 2 ) = A X (N*,Q 2 ) + Ar-(N\Q 2 ), (9) 
Sx(N*,Q 2 ) = S A (A^Q 2 ) + Sr-(iV*,Q 2 ), (10) 

where A^- C -(N*, Q 2 ) and S^~ C -(N*,Q 2 ) are due to the meson cloud effect defined by the 
second term of Eq. PJ . 

With the hadronic parameters of the employed dynamical coupled-channels model de- 
termined in analyzing the irN reaction data 0, @], the only freedom in analyzing the elec- 
tromagnetic meson production reactions is the electromagnetic coupling parameters of the 
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FIG. 2: Fit to p(e,e'7r°)p structure functions at Q 2 = 0.4 (GeV/c) 2 . Here 6 = 9*. The solid curves 
are the results of Fitl, the dashed curves are of Fit2, and the dotted curves are of Fit3. (See text 
for the description of each fit.) The data are taken from Refs. ldl. 111!]. 



model. If the parameters listed in Ref. 0] are used to calculate the non-resonant interaction 
v L'S'M'B',x \ N (k'i O) m Eqs. and (j4j), the only parameters to be determined from the data 
of pion electroproduction reactions are the bare helicity amplitudes denned by Eq. (jSJ) . Such 
a highly constrained analysis was performed in Ref. j3] for pion photoproduction. It was 
found that the available data of jp —>■ n°p, ir + n can be fitted reasonably well up to invariant 
mass W < 1.6 GeV. In this work we extend this effort to analyze the pion electroproduction 
data in the same W region. 

We first try to fix the bare helicity amplitudes by fitting to the data of (Jt + cctl, <jlt, and 
ott of v{ e i e'7r°)p in Ref. [l(| which covers almost all (W, Q 2 ) region we are considering (see 
Table. HJ). In a purely phenomenological approach, we first vary all of the helicity amplitudes 
of 16 bare N* states, considered in analyzing the ttN — > 7tN,7T7tN data 0,0], in the fits 
to the data. It turns out that only the helicity amplitudes of the first iV* states in Sn, 
Pii) -P33 an d -D13 are relevant in the considered W < 1.6 GeV. Thus in this paper only the 
bare helicity amplitudes associated with those four bare N* states (total 10 parameters) are 
varied in the fit and other bare helicity amplitudes are set to zero. The numerical fit is 
performed at each Q 2 independently, using the MINUIT library. 

The results of our fits are the solid curves in the top three rows of Figs. [20 Clearly our 
results from this fit agree with the data well. We obtain similar quality of fits to the data 



of Ref. 10 



at other Q 2 values listed in Table. [I] We have also used the magnetic Ml form 
factor of 7*iV — > A(1232) extracted from previous analyses as data for fitting. The results 
are shown in Fig. [5l We refer the results of this fit to as "Fitl". 

In Fig. [6j we present the G* M , G* E , and G* c form factors of 7*iV — > A(1232) transition 
obtained from Fitl (solid points). In the same figure, we also show the meson cloud effect 
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FIG. 3: Fit to p(e,e'7r°)p structure functions at Q 2 = 0.9 (GeV/c) 2 . Here 9 = 9*. The data are 



taken from Ref. 



Within our model, it has a significant contribution at low Q 2 , but 
2 increases, particularly for G* E and G* c . These results are similar to 

3. 



in the form factors, 
rapidly decreases as 
the previous findings 

The helicity amplitudes of Six, Pn, and D13 resulting from Fitl are shown in Fig. [7J The 
solid circles are the absolute magnitude of the dressed helicity amplitudes (jUJ) and ffTOl . The 
errors there are assigned by MIGRAD in the MINUIT library. More detailed analysis of the 
errors is perhaps needed, but will not be addressed here. The meson cloud effect (dashed 
curves), as defined by A™- c - and S"™' 0, of Eqs. (jHJ) and ( TTOT) and calculated from the second 
term of Eq. (14"]) , are the necessary consequence of the unitarity conditions. They do not 
include the bare helicity term determined here and are already fixed in the photoproduction 
analysis [a]. Within our model (and within Fitl), the meson cloud contribution is relatively 
small in Su and Ai/ 2 of D13 even in the low Q 2 region. 

Here we note that our helicity amplitudes defined in Eqs. ([9]) and (jTUl) are different from 
the commonly used convention, say A c ^ v and Sf™, which are obtained from the imaginary 
part of the 7*iV — > ttN multipole amplitudes [3]. This definition leads to helicity amplitudes 
which are real, while our dressed amplitudes are complex. It was shown in Ref. [l5| that for 
the A(1232) resonance our dressed helicity amplitudes (J2J and (ITU]) can be reduced to A™ 
and 5™ v , if we replace the Green function G^n with its principal value in all loop integrals 
appearing in the calculation. However, such reduction is not so trivial for higher resonance 
states because the unstable ttA, pN,aN channels open, and thus the direct comparison of 
the helicity amplitudes from other analyses becomes unclear. 

At Q 2 = 0.4 (GeV/c) 2 , the data of all structure functions both for p(e,e'ii )p and 
p(e, e'7f + )n are available as seen in Table. [TJ To see the sensitivity of the resulting helic- 
ity amplitudes to the amount of the data included in the fits, we further carry out two fits 
at this Q 2 , referred to as Fit2 and Fit3, respectively. Fit2 (Fit3) further includes the data 
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of Refs. [ll|, ll2j, ll3j (Ref. [11]) in the fit in addition to those of Ref. [lOj which are used in 
Fitl. This means that Fit2 includes all available data both from p(e, eV°)p and p(e, e f ir + )n, 
whereas Fit3 includes the same data but from p(e, eV°)p only. The results of each fit are 
the dashed and dotted curves in Fig. [2]for p(e, e'ir )p and Fig. [S]for p(e, e'ii + )n, respectively. 

The resulting bare helicity amplitudes are listed in the third (Fit2) and fourth (Fit3) 
columns of Table II and compare with that from Fitl. The corresponding change in the 
7iV — > A(1232) form factors and the dressed helicity amplitudes are also shown as open 
circles and triangles in Figs. [6] and [71 A significant change among the three different fits is 
observed in most of the results except G* M in P 33 . This indicates that fitting the data listed 
in Table [I] are far from sufficient to pin down the 7*iV — > N* transition form factors up to 
Q 2 = 1.45 (GeV/c) 2 . It clearly indicates the importance of obtaining data from complete or 
over-complete measurements of most, if not all, of the independent p(e, e'7t)N polarization 
observables. Such measurements were made by Kelly et al. [13] in the A (1232) region and 
will be performed at JLab for wide ranges of W and Q 2 in the next few years [l| . 

It has been seen in Fig. [H]that all of our current fits underestimate <jt of p(e, e'n + )n at 
forward angles. We find that this can be improved by further varying the S31 and Pi 3 bare 
helicity amplitudes within their reasonable range. In Fig. [9j the results with the nonzero S31 
and P13 bare helicity amplitudes (solid curves) are compared with the results without varying 
those amplitudes (dashed curves). The resulting values of the bare helicity amplitudes are 
(A^,Sfy) = (121.6,59.6) and (A^A^Sfy) = (-73.2,-42.9,41.5). The parameters of 
Fit 2 are used for Su, Pn, P 33 , and P i3 in both curves. In the figure we have just shown 
the results at W ~ 1.3 GeV. We confirm that the same consequence is obtained also at 
other W, and find that the P 13 (S31) has contributions mainly at low (high) W. We also 
find that the inclusion of the bare S31 and Pi 3 helicity amplitudes does not change other 
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FIG. 4: Fit to p(e,e'ir°)p structure functions at Q 2 = 1.45 (GeV/c) 2 . Here 6 = 9*. The data are 
taken from Ref. 10]. 
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FIG. 5: G* M normalized by the dipole factor G D = [1 + Q 2 /0.71(GeV/c) 2 ]- 2 . The solid black 
circles at Q 2 > are our fit to the values extracted from previous analyses (those values are taken 
from Ref. 0]). The triangle at Q 2 = is from our photoproduction analysis 0]. 
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FIG. 6: The 7*iV — ► A(1232) form factors. Solid points are from Fitl; dashed curves are the 
meson cloud contribution. Open circles and triangles at Q 2 = 0.4 (GeV/c) 2 are from Fit2 and 
Fit3, respectively. The three points are almost overlapped in G* M . The solid point at Q 2 = is 
obtained from the photoproduction reaction analysis in Ref. 0]. 



structure functions than <jt of p(e, e'iT + )n (at most, most of the change is within the error). 
This indicates that those two helicity amplitudes are rather relevant to p(e, e'ii + )n, but 
not to p(e, e'ir°)p. As shown in Table [U however, no enough data is currently available for 
p(e,e'ir + )n above Q 2 = 0.4 (GeV/c) 2 . The data both of the p(e,e'ii )p and p(e,e'n + )n at 
same Q 2 values are desirable to pin down the Q 2 dependence of the S 3 i and Pi 3 helicity 
amplitudes. 

We now turn to show the coupled- channels effect. In Fig. [TU1 we see that when only 
the tiN intermediate state is kept in the M'B' summation of the non-resonant amplitude 
[Eq. (j5J)] and the dressed 7*iV — > N* vertices [Eq. (j3j)], the predicted total transverse and 
longitudinal cross sections ot and ol of p(e, e'ir )p are changed from the solid to dashed 
curves. This corresponds to only examining the coupled-channels effect on the electromag- 
netic (Q 2 -dependent) part in the 7*iV — ► ttN amplitude. All coupled-channels effects on 
the non-electromagnetic interactions are kept in the calculations. We find that the coupled- 
channels effect tends to decrease when Q 2 increases. This is rather clearly seen in ax- 111 
particular, the coupled-channels effect on o T at high W ~ 1.5 GeV is small (10-20%) already 
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FIG. 7: Extracted helicity amplitudes for S u at W = 1535 MeV (upper panels), P u at W = 1440 
MeV (middle panels), and D\% aXW = 1520 MeV (lower panels). The meaning of each point and 
curve is same as in Fig. [6j 



TABLE II: Ambiguity of resulting bare helicity amplitudes [the results are at Q 2 = 0.4 (GeV/c) 2 ]. 
The errors are assigned by MIGRAD in the MINUIT library. 
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100.80 ± 1.46 


83.25 ± 1.21 
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-119.30 ±20.41 


-9.85 ± 1.69 


-53.53 ±4.75 


Pu A 1/2 


33.18 ±2.11 


-15.68 ±1.00 


20.17 ± 10.37 


Pll S\/2 


37.29 ± 2.26 


52.23 ±3.16 


131.00 ±5.87 


-P33 ^3/2 


-146.00 ±0.60 


-137.50 ±0.56 


-150.80 ± 1.03 


^33 A1/2 


-54.47 ±0.61 


-62.57 ±0.69 


-46.29 ± 1.73 


-P33 S\/ 2 


7.85 ± 1.25 


-7.66 ±1.22 


7.34 ± 1.69 




-44.01 ± 1.31 


-67.01 ± 1.99 


-98.63 ±2.92 


D13 A 1/2 


97.11 ±8.51 


14.34 ±1.26 


70.02 ± 4.83 


Dl3 5*1/2 


-18.35 ± 1.37 


19.43 ± 1.45 


4.11 ±2.76 



at Q 2 = 0.4 (GeV/c) 2 . (The effect is about 30-40% at Q 2 = [4].) This is understood as 
follows. In Eq. we can further split the resonant amplitude t R as t R = + where 
t R are an d ^m.c. are the same as t R but replacing f J N * A Ajv with its bare part T J N , X Ajy and meson 
cloud part [the second term of Eq. (HJ)], respectively. The coupled-channels effect shown in 
Fig. [10] comes from tj jSNnNX ^ XN and t R c - We have found that the relative importance of the 
coupled-channels effect in each part remains the same for increasing Q 2 . However, the contri- 
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FIG. 8: Structure functions of p(e, e'ir + )n at Q 2 = 0.4 (GeV/c) 2 . Here 6 = 9*. The solid curves 
are the results of Fitl, the dashed curves are of Fit2, and the dotted curves are of Fit3. (See text 
for the description of each fit.) As for the glt'i results at W = 1.14, 1.22, 1.3, 1.38, 1.5, 1.58 GeV 
(from left to right of the bottom row) are shown, in which the data are available. The data in the 
figure are taken from Ref . 12J, [13 ] . 



bution of non-resonant mechanisms both on t J Ls N ^N,Xy\ N an< ^ *m.c. to the structure functions 
decreases for higher Q 2 compared with t^ am . This explains the smaller coupled-channels 
effect compared with the photoproduction reactions [3]. The decreasing non- resonant inter- 
action at higher Q 2 is due to its long range nature, thus indicating that higher Q 2 reactions 
provide a clearer probe of N*. We obtain similar results also for p(e, e'7t + )n. 

It is noted, however, that the above argument does not mean coupled-channels effect 
is negligible in the full 7*iV — > ttN reaction process. In the above analysis we kept the 
coupled-channels effect on the hadronic non-resonant amplitudes, the strong N* vertices, 
and the N* self-energy, which are Q 2 -independent and remain important irrespective of Q 2 . 
We have found in the previous analyses [1, 0| that the coupled-channels effect on them is 
significant in all energy region up to W = 2 GeV. 

In Fig. [TTJ we show the coupled-channels effect on the five-fold differential cross section 
defined by Eq. (JSJ). The coupled-channels effects are significant at low W, whereas they 
are small at high W . This is consistent with the above discussions because the five-fold 
differential cross sections are dominated by cry. Here we also see that our full results (solid 
curves) are in good agreement with the original data, although we performed the fits by 
using the structure function data listed in Table [D 
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FIG. 9: Contribution of the S31 and P 13 helicity amplitudes at Q 2 = 0.4 (GeV/c) 2 . The left (right) 
panels are the structure functions of p(e, e'7r°)p \p(e,e'ir + )n] reaction at W = 1.3 GeV (W = 1.29 
GeV). Solid (dashed) curves are the results with (without) nonzero and P13 bare helicity 
amplitudes. The parameters of Fit2 are used for the Su, Pn, P33, and D13 helicity amplitudes in 



both curves. The data are from Refs. 
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III. SUMMARY AND OUTLOOK 

In this work we have explored how the available p(e, e'ir)N data can be used to determine 
the 7*A" — > N* transition form factors within a dynamical coupled- channels models 0, 0, 
0, 0, 0] • Within the available data, the 7*iV — > A^* bare helicity amplitudes of the first N* 
states in ^n, Pn, P33 and Di 3 can be determined in the considered energy region W < 
1.6 GeV. We further observe that some of these parameters can not be determined well. 
The uncertainties could be due to the limitation that only data of 4 out of 11 independent 
p(e, e'n)N observables are available for our analysis. Clearly, the data from the forthcoming 
measurements of double and triple polarization observables at JLab will be highly desirable 
to make progress. 

Also, it was found that the underestimation of the <jt of p(e, e'ir + )n at forward angles can 
be improved by further considering the S31 and P13 bare helicity amplitudes. Furthermore, 
these amplitudes can have relevant contribution to p(e, e'ir + )n, but not to p(e, eV )p. The 
p(e, e'ir + )n data of wide Q 2 region as well as p(e, eV°)p seem necessary for determining the 
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FIG. 10: Coupled-channels effect on the integrated structure functions gt{W) and <jl{W) for Q 2 = 
0.4, 0.9, 1.45 (GeV/c) 2 for p(e,e'TT°)p reactions. The solid curves are the full results calculated 
with the bare helicity amplitudes of Fitl. The dashed curves are the same as solid curves but only 
the irN loop is taken in the M'B' summation in Eqs. (|2|) and (JH). 



Q 2 dependence of the S31 and P13 helicity amplitudes. 

For testing theoretical predictions from hadron structure calculations such as LQCD, the 
quantities of interest are the residues of the 7*iV — ► nN amplitudes, denned by Eqs. ((!])- (jlj), 
at the corresponding resonance poles. If the resonance poles are associated with the am- 
plitude t^g N7rNX \ N (k, q, W, Q 2 ) of Eq. (j3J), the extracted residues are directly related to the 
dressed form factors T J N * l'S'M'B' (^'1 Wl: An analytic continuation method for extracting 
these information has been developed [181 ]. and our results along with other hadronic prop- 
erties associated nucleon resonances will be published elsewhere. Here we only mention that 
the extracted form factors are complex and some investigations are needed to see how they 



can be compared with the helicity amplitudes, which are real numbers, listed by PDG [19 
In a Hamiltonian formulation as taken in our dynamical approach, the physical meanings of 
poles and residues are well defined in textbooks (HHj. 
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FIG. 11: Coupled-channels effect on the five- fold differential cross sections T^lda 5 /(dE^dn^dn^ )] 
of p(e,e'ir°)p (upper panels) and p(e, e'ir + )n (lower panels) at Q 2 = 0.4 (GeV/c) 2 . Here 9 = 9* 
and (ft = (ft*. The solid curves are the full results calculated with the bare helicity amplitudes of 
Fitl. The dashed curves are the same as the solid curves but only the ttN loop is taken in the 



M'B' summation in Eqs. ([2]) and (j3|). The data are taken from Ref. 14] 



DE-AC02-05CH11231. 
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